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Abstract — To remove random valued impulse noise from 
colour images, an efficient impulse detection and filtering 
scheme is presented. The locally adaptive threshold for 
impulse detection is derived from the pixels of the filtering 
window. The restoration of the noisy pixel is done on the basis 
of brightness and chromaticity information obtained from the 
neighbouring pixels in the filtering window. Experimental 
results demonstrate that the proposed scheme yields much 
superior performance in comparison with other colour image 
filtering methods. 

Index Terms — Impulse detector, image filtering, random 
valued impulse noise. 

I. Introduction 

Image filtering is an essential part of any image 
processing system. A large number of filtering techniques 
have been developed for removal of impulse noise from 
colour images [1-2]. Most of these techniques use vector 
processing approach as it is widely accepted that this 
approach is more appropriate than the component-wise 
filtering approach, which can generate colour artefacts in 
the filtered image. The vector median filter (VMF) [3], 
vector directional filter (VDF) [4] and the directional 
distance filter (DDF) [5] are the most commonly used 
filters for noise removal in colour images. However, if the 
filtering operation is preceded by an impulse detection 
stage separately, the component-wise approach for impulse 
detection can be used with good results. The impulse 
detection stage can be followed by the filtering operation 
using the vector approach. The presence of an impulse 
detector ensures that only the noisy pixels are filterd. This 
prevents the blurring which is caused by the filtering of 
noise-free pixels. Indeed, some of the recent methods such 
as high performance detection (HPD) [6], directional 
weighted median filter (DWM) [7], and switching vector 
median filter based on CIELAB colour space [8] are based 
on separate impulse detection scheme. Filtering of only 
noisy pixels results into better preservation of detail 
features. 

The main drawback of applying the component-wise 
approach for filtering is that the inherent correlation 
existing among pixels of different channels may be lost, 
resulting into colour artefacts. To overcome this problem, 
color of filtered pixel may be corrected separately. In this 
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way, the ability of component-wise filtering approach to 
produce lower mean square error can still be utilized if it is 
combined with some color correction scheme. 

In this paper we present an impulse detection method 
where threshold is adaptively changed on the basis of local 
pixel statistics within the filtering window. This is followed 
by a filtering method which takes into consideration the 
colour and brightness information to restore noisy pixels. 
The proposed method yields much better peak signal to 
noise ratio (PSNR), mean absolute error(MAE) and 
normalized colour difference (NCD) than other methods as 
it combines the best features of both component-wise 
processing and vector processing approaches. 

II. Impulse Detection 

The impulse detection is based on the assumption that 
a noisy pixel takes a value which is substantially different 
than the neighbouring pixels in the filtering window, 
whereas noise-free regions in the image have locally 
smoothly varying values separated by edges. 

LetXj . = ixf], xfj , xf' J be a multi-channel pixel in the 
RGB space at location (i, j) of image /, which is corrupted 
by random valued impulse noise. Let us also define yi/'J, 
as (w x n) filtering window for channel k with centre pixel 






\ . Impulse detection is performed for each channel 



separately. First of all, the median, med (w ( 5 k J 5 ) , is 
subtracted from each pixel in the window w^] to obtain 
the differences as: 

w%ff=w™-med(yv£! i ) ; k =1,2,3 (1) 

Now for each value k, the differences in the window 
wy ff are arranged in ascending order as 

\dL', d\J, , d, 2 L | . Then a parameter r t) , defining the 

roughness of the filtering window is computed as: 



'« 



"M" 



(2) 



Now an adaptive threshold, which depends on statistical 
characteristics of pixels within the window w**, for each 
channel, is empirically obtained for natural images as: 

T (k) =15 + 2.6r w xexp(-0.003r w 2 ) (3) 
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The output of the detector is represented in terms of a 
flag image If, , 1 where 

\l ; if 4 k j-med(wtl)>T {k) 




Here / ; y = 1 implies that x\ is noisy. 



A. Image Restoration 

In the first phase of filtering to restore the brightness, a 
directional weighted scheme, as used in [7], is considered. 
The weight of a pixel is decided on the basis of standard 
deviation in four pixel directions (vertical, horizontal and 

two diagonals). The brightness of the noisy pixel I f}y = 1) 
is restored as: 



yfl =med\w lm 0x\ 



.(*) 



for all / and m such that 



where operator denotes repetition operation and weights 
{ w, m } are defined as: 

\2-MUs (5) 

[l ; otherwise 

Here S denotes the set of pixels in the direction with 
minimum standard deviation. Now, the brightness restored 

multi-channel pixel is represented as u, , = I up, , uf', , u\\ I 

where 



.(*) 



.,(*) 



(*) 



tfytf/tf=i 



M. 



otherwise 



for k =1,2,3 



(6) 



It is well known that VDF can provide the best possible 
color information about the filtered pixel. However, the 
magnitude of the filtered pixel vector is not very reliable in 
this case. Therefore, we use the colour information 
provided by the VDF in the form of direction of the filtered 
pixel vector, along with the magnitude information 
obtained by (6) for better results. In the second phase of 
filtering, colour of the x, . is restored by using the colour 

information obtained from VDF. The output 

V. . =| v' 1 ', vr'.v' 3 ' 1 of the VDF is obtained from a 

3x3filtering window by considering noise free-pixels. 
A 5 x 5 is window is used if the number of noise-free pixels 
is less than four in 3x3 window. The final output after 

colour restoration is given as 2^=1 z\ ) , 

f (*) II II i\\ II -c A*) 1 
' 'J 'J ' fl.J 



,Zi!,Zij) where 



#~ 



,« 



otherwise 



(7) 



correlated random valued impulse noise with noise density 
varying from 5% to 30%. The criteria used to compare the 
performance of various filters are PSNR (dB), MAE and 
NCD. These criteria are defined as: 



PSNR = 1 log, 



(1/3MA02; I lh.j-o 

M N 

MAE = (l/3MW)^^|z ij -o,. 



(8) 



(9) 



where M and N define the image size, zy is the pixel from 
the filtered image and Oy is the corresponding pixel in the 

original image. The term ||.||~ denotes the Euclidian distance 
and |.| denotes the City-block distance. For NCD 

calculations, first of all, RGB to CIELAB conversion 
[1,2,9] is done and then we use the following expression: 
NCD = 

M N / i , 

Ill[£(u)-n(iJ)] + [<(y)-«;(u)] +[*;(uH:(u)J 

M N j j 9 7 \l/2 

Ili[n(u)] + [<c(i,j)] + [K(ij)]) 

(10) 
where a* (i, j") , b* (/, j) denote the chrominance value and 

L* (;', y) denotes the lightness value of the filtered pixel in 
the CIELAB colour space. These parameters for noise free 
pixel of the original image are a* (i, _/') ,b* (i, j) and 

L* ((', j) , respectively. 

The restoration results of various filters for subjective 
visual qualities are shown in the Fig. 1 for 'Peppers' image. 
The PSNR, MAE and NCD resulting from various 
experiments are presented in Table I and II for 'Peppers' 
and table III and IV 'Boat in Lake' images, respectively. 
From these tables, it is observed that the proposed method 
outperforms the other methods at all noise levels in terms 
of PSNR, MAE and NCD. The colour correction of noisy 
pixels results in the lowest NCD values as evident from the 
results. The self adjusting nature of the threshold ensures 
the reliable detection of noisy pixels. The detection is 
followed by an equally effective filtering method that 
corrects both brightness and colour for restoration of noisy 
images. 



III. Simulation Results 

To assess the performance of the proposed method, we 
compare it with several methods including VMF, VDF, 
DWM, HPD and that of Jin et al. [8]. The test images used 
in simulations are 'Peppers' and 'Boat in Lake' each of size 
512x512. The test images are corrupted with 50% 
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TABLE I. 

Comparison of restoration results in PSNR(db) for 
"Peppers' image 



Method 


Noise percentage 


5 


10 


15 


20 


25 


30 


VMF 


35.8 


34.4 


32.7 


30.1 


27.6 


25.4 


VDF 


34.4 


32.7 


30.5 


28.0 


25.6 


23.2 


Jin [8] 


38.7 


35.8 


33.4 


30.6 


28.0 


25.5 


HPD 


38.1 


36.5 


35.2 


34.1 


33.4 


31.0 


DWM 


38.9 


37.5 


36.1 


34.3 


32.2 


31.1 


Proposed 


40.5 


38.7 


37.0 


34.7 


33.9 


31.3 



TABLE II. 
Comparison of restoration results in mae & ncd for 
'peppers' image 



Method 


Noise percentage 


10 


20 


MAE 


NCD 


MAE 


NCD 


VMF 


2.22 


.0021 


3.29 


.0080 


VDF 


2.63 


.0022 


4.06 


.0091 


Jin [8] 


0.98 


.0017 


2.25 


.0077 


HPD 


0.61 


.0019 


1.12 


.0036 


DWM 


0.56 


.0013 


1.05 


.0027 


Proposed 


0.46 


.0012 


1.04 


.0026 



TABLE III. 

comparison of restoration results in psnr(db) for 'boat in 
Lake' image 



Method 


Noise percentage 


5 


10 


15 


20 


25 


30 


VMF 


30.9 


30.0 


29.0 


27.5 


25.8 


23.8 


VDF 


29.1 


28.2 


27.2 


25.6 


24.5 


22.7 


Jin [8] 


34.9 


32.6 


30.5 


28.9 


26.4 


24.3 


HPD 


32.6 


31.4 


30.3 


29.5 


28.7 


28.0 


DWM 


34.8 


32.7 


31.1 


30.2 


29.3 


27.9 


Proposed 


35.2 


33.6 


32.9 


31.2 


29.6 


28.2 



TABLE IV. 
Comparison of restoration results in mae & ncd for 'Boat 
in Lake' image 



Method 


Noise percentage 


10 


20 


MAE 


NCD 


MAE 


NCD 


VMF 


4.01 


.0025 


5.41 


.0073 


VDF 


4.06 


.0026 


5.28 


.0067 


Jin [8] 


1.67 


.0018 


3.51 


.0063 


HPD 


1.45 


.0038 


2.31 


.0068 


DWM 


1.66 


.0089 


.28 


.0178 


Proposed 


0.99 


.0017 


1.79 


.0042 



Conclusions 

We propose an impulse detection and filtering scheme to 
remove random valued impulse noise from colour images. 
The detection of noisy pixels is based on a adaptive 
threshold which is dependent on the local statistics in the 
filtering window. The filtering method corrects only the 
corrupted pixel values by taking into consideration the 
brightness and colour information for restoration of noisy 
images. The efficacy of the proposed method is 
demonstrated by experimental results, which exhibit 
significant improvement over several other methods. 




Figure 2. Restoration performance of different filters for Peppers 

image. (a) Original (b) Noisy with noise density 25% (c) VMF (d) VDF 

(e) Jin et al[6] (f) DWM (g) HPD (h) Proposed. 



References 

[1] R C Gonzales and R E Woods, Digital Image Processing, 2" 

ed. Reading, MA: Addison -Wesley, 2002. 
[2] R. Lukac and K. N. Plataniotis, Color Image Processing: 

Methods and Applications , Taylor and Francis, 2007. 
[3] J. Astola, P. Haavisto and Y.Neuov, "Vector median filters", 

Proc. IEEE, 1990, 78, (4) pp.678-689. 
[4] P. E. Trahanias, and A. N. Venetsanopoulos, "Vector 

directional filters - a new class of multichannel image 

processing filters", IEEE Trans. Image Process., 1993, 2, ( 

4), pp. 528 - 534. 
[5] D. G. Karakos, and P. E. Trahanias, "Generalized 

multichannel image-filtering structures", IEEE Trans. Image 

Process., 1997, 6, (7), pp. 1038 - 1045. 
[6] A.S. Awad and H. Man, "High performance detection filter 

for impulse noise removal in images", Electronics Letters, 

vol. 44, issue 3 pp.192-194 , Jan 2008. 
[7] Y.Dong and S. Xu, "A new directional weighted median 

filter for removal of random-valued impulse noise", IEEE 

Signal Processing Letters, vol. 14, issue 3, pp. 193-196, 

March 2007. 
[8] L. Jin and D.Li, "A switching vector median filter based on 

the CIELAB color space for color image restoration", Signal 

Process., 2007, 87, (6), pp. 1345-1354 
[9] K. Jack, Video Demystified: A Handbook For The Digital 

Engineer, Elsevier, 2007. 



©2010 ACEEE 

DOI: 01.IJSIP.01.03 .235 



;C 



ACEEE 



